% Marr-Hildreth (LoG) 检测
sigma = 2; % 高斯模糊标准差
hsize = 2 * ceil(3 * sigma) + 1; % 窗口大小
log_filter = fspecial('log', hsize, sigma); % 生成 LoG 滤波器

% 应用 LoG 滤波
log_response = imfilter(double(enhanced_image), log_filter, 'replicate');

% 阈值化得到边缘
log_edges = log_response < -0.1; % 根据实验调整阈值
imwrite(log_edges, 'log_edges.png');
